Comparative study of the electronic structures of Fe 3 04 and Fe 2 Si04 
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The electronic properties of two spinels FesCU and Fe2SiC>4 are studied by the density functional 
theory. The local Coulomb repulsion U and the Hund's exchange J between the 3d electrons on 
iron are included. For U — 0, both spinels are half-metals, with the minority ti g states at the Fermi 
level. Magnetite remains a metal in a cubic phase even at large values of U. The metal-insulator 
transition is induced by the X3 phonon, which lowers the total energy and stabilizes the charge- 



orbital ordering. Fe2SiC>4 transforms to a Mott insulating state for U > 2 eV with a gap A 9 



U. 



The antiferromagnetic interactions induce the tetragonal distortion, which releases the geometrical 
frustration and stabilizes the long-range order. The differences of electronic structures in the high- 
symmetry cubic phases and the distorted low-symmetry phases of both spinels are discussed. 



I. INTRODUCTION 

In transition-metal spinels, magnetic, charge, and or- 
bital interactions are strongly affected by the crystal ge- 
ometry. A special character of spinels is associated with 
the arrangement of cations in the octahedrally coordi- 
nated sites forming a pyrochlore lattice. Consisting of the 
corner-sharing tetrahedra, this crystal structure causes 
the geometrical frustration of magnetic interactions and 
makes the electronic ground state highly degenerate [lj. 
In such cases, the system may remain in a disordered 
state like a quantum spin liquid or spin ice even at low 
temperatures [2|- It makes also the electronic state very 
susceptible to lattice distortion, which lifts the degener- 
acy and drives the structural phase transition. 

Spinels have a general formula AB2O4 and at room 
temperature crystallize in the cubic Fd3m (O^) symme- 
try (Fig. [T|). In a normal spinel, the A and B cations 
occupy the tetrahedral and octahedral positions, respec- 
tively. The B sites create the atomic chains along the 
[110] and [110] directions. In the most common cases, 
the valence state for the A site is 2+ and for the B site 
3+ (e.g. in MgA^C^). The 4+ valence state for the 
A site exists only in the high-pressure forms of olivines 
(Fe2Si04, Co 2 Si04). In inverse spinels, the A atoms oc- 
cupy half of the octahedral sites, while the B ones are 
shared equally by the octahedral and tetrahedral sites. 
Fe3C>4 is classified as the inverse spinel due to mixed va- 
lency of iron ions in the octahedral sites (2+/3+). 

Depending on the type of cations in the octahedral 
sites, transition-metal spinels exhibit various types of or- 
dered phases at low temperatures [3| . One observes anti- 
ferromagnetic (AF) phases (MgV^CM) [J], charge order- 
ing (CUIX2S4) Q, metal-insulator transitions (MgTi2C>4) 
@, or spin-Peierls transitions (ZnCr 2 04) 0. Crystal dis- 
tortions in these phase transitions are associated with dif- 
ferent electron- lattice couplings. In ZnMn204, the sym- 
metry reduction to the tetragonal space group IAi/amd 
is driven by the Jahn- Teller effect. The tetragonal dis- 
tortion in ZnCr2C>4 observed at Tjv — 12.5 K induces the 
spin gap in magnetic excitations similarly as in the one- 
dimensional spin-Peierls systems. The metal-insulator 
transition in MgTi2C>4 is associated with the spin-singlet 



formation and dimerization of the atomic bonds [8( . 

The Verwey transition in Fe3C>4 has a very unique 
character and its microscopic origin has been intensively 
studied for more than 60 years. Verwey proposed that 
a discontinuous drop of the electrical conductivity at 
Ty = 122 K is associated with the ordering of the Fe 3+ 
and Fe 2+ ions in the octahedral sites Q. According to 
this scenario, the conductivity increases above Ty due 
to higher mobility of electrons fluctuating between the 
Fe sites. At Ty, electron movement freezes out and the 
charge ordering stabilized by the electrostatic forces de- 
velops. However, this simple picture was not confirmed 
by more recent studies. Firstly, Anderson demonstrated 
that in the spinel geometry the energy difference be- 
tween a short-range and long-range order is very small 
and therefore typical order-disorder mechanism is there- 
fore unlikely [lj. Secondly, the crystallographic studies 
revealed the monoclinic deformation below Ty, which is 
inconsistent with the Verwey model [lCj] • The diffraction 
measurements support only the fractional charge order- 
ing with small differences of charges [TJj UM ■ Finally, the 
theoretical studies demonstrated the importance of local 
Coulomb interactions in the 3d states on iron, which in- 
duces the charge-orbital ordering in the monoclinic phase 




FIG. 1: The spinel structure: the octahedral B sites are rep- 
resented by red balls and the AO4 units are plotted as green 
tetrahedra. 



|13j . In our studies, we developed the mechanism of the 
Verwey transition based on the interplay between the 
electron correlations and the electron-phonon interaction 

dHH. 

Spinels containg Si atoms belong to a separate group 
of minerals stable only under high pressure, as real- 
ized in the Earth upper mantle. Fe2Si04, which is 
the end-member of geologically import group of olivines 
(Mg,Fe)2SiC>4, can be obtained from magnetite by re- 
placing the Fe atom at the A site by the Si atom. The 
experimental information on the electronic properties of 
this mineral is very limited. The heat capacity measure- 
ment revealed the peak at T/v = 11.8 K, most proba- 
bly associated with the paramagnetic-antiferromagnetic 
transition [17J. Above Xjv, a paramagnetic state is consis- 
tent with the Mossbauer measurements [18(. In the pre- 
vious study, we have analyzed the electronic structure of 
Fe2SiC>4 using the density functional theory (DFT) [19j . 
We found that the insulating state results from the local 
Hubbard interaction U included within the generalized 
gradient approximation within the so-called GGA+U ap- 
proach. In addition, we revealed that the AF interactions 
induces the tetragonal distortion, removing the geometri- 
cal frustration. In this paper, we compare the electronic 
properties of Fe3C>4 and Fe2SiC>4, analyzing in details the 
changes induced by the Hubbard interaction U and crys- 
tal distortions. For Fe2SiC>4, we investigate the effect of 
high pressure on the electronic structure. 

The paper is organized as follows. The details of the 
calculation method are presented in Sec. HU Electronic 
structures of Fe3C>4 and Fe2SiC>4 are reported and ana- 
lyzed separately for the high symmetry phase (Sec. IIII Ap 
and low symmetry phase (Sec. IIII B[) . Discussion of the 
results and a short summary are presented in IIVI 



II. CALCULATION METHOD 

The electronic and crystal structures of both spinels 
were optimized within the GGA+U method using the 
VASP code [20]. The wave functions were obtaine d by 
the full-potential projector augmented- wave method |2l| . 
The basis included the following valence electron configu- 
rations: Fe:3d 6 4s 2 , Si:3s 2 3p 2 , and 0:2s 2 2p 4 . To compare 
the effects induced by the local electron interactions in 
the 3d states on iron, the same values of the Coulomb and 
Hund's exchange parameters were chosen: U — 4 eV and 
J = 0.8 eV. The repulsion energy U was taken from the 
constrained DFT calculations 23] and J was obtained 
from the atomic values of the Racah parameters for Fe 2+ 
ions [H: B = 0.131 eV and C = 0.484 eV. It is reason- 
able to assume that the Coulomb interaction parameters 
{U, J} are the same in both Fe3 04 and Fe2Si04 systems 
in their insulating phases, as one expects that the screen- 
ing of the Coulomb interaction should be rather similar. 
However, we cannot exclude that the Coulomb parame- 
ter U changes in the magnetite at the Verwey transition, 
when the metallic state is transformed into an insulating 



one (this effect was not considered in the theory until 
now). 

The calculations for the cubic symmetry FdSm were 
performed in the a x a x a supercell and for the mono- 
clinic phase P2/c of magnetite in the o/y2 x a/y/2 x 2a 
cell, both with 56 atoms. The electronic structures were 
obtained for the fully relaxed crystals. For Fe 2 Si04, the 
results were obtained at two hydrostatic pressures p = 
and p = 20 GPa. The latter value corresponds to a stable 
spinel phase existing only under high-pressure conditions. 



III. ELECTRONIC STRUCTURES 
A. High-symmetry phases 

In the case of the iron-based spinels, the electronic 
structures can be well understood starting from the 
electron configurations of the Fe 2+ and Fe 3+ ions. In 
Fe2Si04, the iron atoms have a uniform electron distribu- 
tion corresponding to the 2+ charge valency. The octahe- 
dral sites in magnetite have the mixed valence configura- 
tion with an average charge 2.5+. Due to the Hunds's ex- 
change both ionic configurations correspond to the high 
spin states Fe 2+ : t 3 t t 2g 4, e 2 t (S=2) and Fe 3+ : t 3 , +e 
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(S=5/2). The splitting of the 3d states into the tig and 
e g orbitals is induced by the crystal field, which is much 
larger at the octahedral sites than at the tetrahedral ones. 
The exchange coupling causes the splitting of the major- 
ity up-spin and minority down-spin states and shifts one 
type of the states with respect to the other. Consequently 
only the down-spin states occupy the energy levels close 
the Fermi energy (Ep) j23| . 

In Fig. [5J we compare the electronic density of states 
(DOS) in these compounds for U = J = 0. If the local 
electron interactions are neglected, the grounds states of 
both spinels are half-metallic with the minority down- 
spin t<i g states at Ep. These states come only from the 
Fe atoms at the octahedral sites. Neither the 3d states 
on the Fe(A) atoms in magnetite nor the s and p states 
participate in the electron transport. Similarly to the 
Fe(B) states, the up-spin and down-spin states at the 
Fe(A) sites are splitted by the exchange coupling. All the 
down-spin states are occupied, while the up-spin states 
are empty, so the magnetic moments on the A and B 
sites have the opposite directions generating the ferri- 
magnetic order in magnetite. The half-metallic ground 
state of magnetite agrees with the properties of the high- 
symmetry phase indicating that the local Coulomb inter- 
actions do not play important role at high temperatures. 

Without the Hubbard interaction, the metallic ground 
state of Fe2Si04 has the ferromagnetic (FM) order with 
the magnetic moment on iron equal 3.57/is. It remains 
a metal even at high pressures, showing a disagreement 
with the insulating behavior of this material. As shown 
in Fig. 2(b), the main effect of pressure is to increase the 
bandwidth due to larger overlap of the 3d orbitals. The 
splitting of the t2 g and e g states also increases because 
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FIG. 2: Electronic structure of (a) Fe304 at p = pressure 
and (b) Fe2Si04 obtained in the cubic Fd3m symmetry for 
two pressure values, p = (solid lines) and p = 26 GPa 
(dashed lines), and for U — J — 0. 



of the stronger crystal field. 

In the next step, we compare the changes in the DOSs 
induced by the on-site interactions U and J, while keep- 
ing the cubic symmetry. We allow only for the relax- 
ation of the lattice constant and the internal parameter 
u. The lattice constants of both crystals increase slightly 
(less than 1%) because of the larger Coulomb repulsion 
between electrons [15|,[lj|. ^ n magnetite, the electronic 
state does not change much comparing to the uncorre- 
cted case, Fig. 3(a). There is an additional shift of 
the majority up-spin Fe(B) and down-spin Fe(A) states 
to lower energies due to larger exchange coupling. It 
increases the magnitudes of magnetic moments on the 
Fe(A) sites from 3.40 to 3.97 /ib and on the Fe(B) sites 
from 3.52 and 3.87 /is, improving the agreement with 
the experiment. The ground state remains metallic even 
at presence of strong on-site interaction U . 

In Fe2Si04, there is a larger electron occupation of 
the t2 g states on the Fe(B) atoms (one electron per one 
iron site), so the on-site Hubbard interaction U is much 
more effective than in magnetite. As shown in Fig. 3(b), 
the interaction U has a strong impact on the electronic 
DOS, inducing the insulating gap at Ep. The gap with 
the magnitude 1.5 eV separates the occupied down-spin 
ti g states from the empty ones. The insulating state is 
generated by the inter-orbital Hubbard interaction with 
the effective energy U e ff = U — J (note that U is here an 
average value of local Coulomb interaction for different 



FIG. 3: The same as in Fig. 1 for U = 4 eV and J = 0.8 eV. 



electron pairs, as usually assumed in the local density 
approximation supplemented by local electron-electron 
interactions (LDA+J7) [13|, [24]]), so the magnitute of the 
gap can be estimated from the formula 



A„ 



U, 



eff 



w 

Y' 



(i) 



where W is the width of the ti g band. For W ~ 3 eV, 
we get A g ~ 1.7 eV, which corresponds very well with 
the obtained gap. At high pressure, the insulating gap 
decreases to 1.2 eV, and this effect can be easily under- 
stood from Eq. (JTJ) — under increasing pressure the band 
width W increases, thus reducing the value of the gap. 

The results obtained for a cubic phase with a realis- 
tic value of U reveal a fundamental difference between 
these two spinel structures. On the one hand, Fe2Si04 
behaves like a Mott insulator with the gap A g ~ U, and 
the local Hubbard interaction is crucial to reproduce the 
electronic properties of this material. On the other hand, 
the on-site interaction U alone does not induce the insu- 
lating state in magnetite, and in the cubic symmetry this 
material remains in a metallic state. This may be con- 
cluded from the behavior at high temperatures, in partic- 
ular from a relatively high conductivity. This theoretical 
result indicates that Fe304 is not a Mott insulator, and 
its non-metallic state at low-temperatures has a different 
origin. 



B. Low-symmetry phases 

The phase transitions observed in the discussed ma- 
terials at low temperatures have fundamentally different 
nature. The Verwey transition involves both the charge 
and orbital degrees of freedom, and their ordering below 
Ty is stabilized by the monoclinic distortion 13]. Us- 
ing the group theory, we found that the low-symmetry 
phase results from the condensation of two phonons (or- 
der parameters) with the symmetries A5 and X3 [l4|. 
The strength of the electron-phonon coupling strongly 
depends on the Hubbard interaction U, which induces the 
antiferro-orbital ordering and splitting of the t 2g states 
at the Fermi level. Without the lattice distortion, the 
electronic state with orbital polarization is highly degen- 
erate showing a short-range order [25|. We found that 
the charge-orbital ordering is stabilized by the coupling 
between the 3d electrons and the X3 phonon [lq . This 
coupling lowers the total energy of the crystal and in- 
duces the metal-insulator transition. 

The orbital order of the t^g states with the atom dis- 
placements corresponding to the X3 mode is shown in 
Fig. 4(a). The positions of atoms were generated from 
the polarization vectors using the PHONON software 
[26|. Iron atoms occupy 4 different positions in the stru- 
cure labeled as Bi, B2, B3, and B4. Phonons at the X 
point are doubly degenerate, so we have chosen only one 
phonon branch. In the presented case, the displacements 
of the B3 and B4 atoms together with oxygens located in 
the same plane induce the orbital ordering on the Bi and 
B2 atoms in the nearest neighbor planes. The second or- 
der parameter A5 also couples strongly to electrons but it 
does not induce the gap opening. It generates the crystal 
distortion which doubles the unit cell in the c direction. 
These results demonstrate that the cubic phase is not a 
real ground state of magnetite, and the structural transi- 
tion to a phase with broken symmetry must occur when 
temperature is lowered. 
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FIG. 4: Orbital ordering induced by (a) the X3 mode 
(the black arrows), and (b) the monoclinic symmetry P2/c. 
Nonequivalent Fe positions are labeled BI, B2, B3, and B4. 




SiO. 



FIG. 5: The AF ordering in Fe 2 Si04 — light/dark spheres 
indicate positions of t and J, magnetic moments on Fe atoms. 



In contrast, the phase transition occurs in Fe2Si04 at 
much lower temperature TV = 11.8 K, and is driven 
by the magnetic interactions. We have considered two 
magnetic configurations with the FM and AF order. We 
assumed that the magnetic moments in the AF configu- 
ration are aligned parallely in the (001) planes, and have 
opposite (alterneting) directions in the subsequent planes 
along the c axis (see Fig. [5]). Instead of the c axis one 
can choose either a or b direction, so the AF order is 
triply degenerate in the cubic symmetry. We found that 
for U > 2 eV, there is a change of the ground state from 
the FM to AF order [lj|. This change is driven by the 
metal-insulator transition at the same value of U. There- 
fore, the electronic DOS with the FM order presented in 
Fig. 3(b) does not correspond to the ground state real- 
ized at U = 4 eV. We revealed that the AF configuration 
induces the tetragonal distortion (a > c), breaking the 
cubic symmetry and lowering the total energy. This dis- 
tortion depends on the Hubbard interaction U and for 
U = 4 eV one finds c = 0.985a. The contraction of the 
crystal along the c axis enhances the superexchange in- 
teraction along the AF bonds, while the elongation in the 
ab plane weakens the one along FM bonds. 

The electronic DOSs in the low symmetry phases of 
both spinels were obtained by full relaxation of the crys- 
tal structures without imposing the cubic symmetry con- 
strains. Fe304 has the monoclinic P2/c symmetry with 
the doubled lattice constant in the c direction. In this 
symmetry there are two non-equivalent tetrahedral posi- 
tions Ai and Ai (not shown), and four octahedral sites 
Bi, B 2 , B 3 , and B4, shown in Fig. 4(b). The lattice 
parameters and atomic positions were analyzed in detail 
in Ref. [15(. Fe2Si04 has the tetragonal Mi/amd space 
group. 

The electronic structures for the above symmetries are 
presented in Fig. \5\ The insulating gaps in the minority 
t2 g states now exist in both systems. Comparing to the 
FM phase of Fe2Si04, the gap in the AF phase is larger (2 
eV). Under pressure, the magnitude of the gap decreases 
as in the FM case. The gap A g = 0.33 eV obtained in 
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FIG. 6: Electronic structures of (a) the P2/c phase of FesCU, 
and (b) the IAi/amd phase of Fe2Si04. Parameters: U = 4 
eV and J = 0.8 eV. 



magnetite agrees quite well with the experimental value 
0.14 eV [27| and is much smaller than the gap in Fe2Si04. 
The gap separates the occupied states at the Bi and B4 
sites from the empty ones at the B2 and B3 sites. This 
explains the origin of the fractional charge ordering in the 
monoclinic phase. The orbital ordering arises due to the 
preferential occupation of the tig states. A schematic plot 
of the orbital ordering is presented in Fig. 4(b). In the 
Bi chains, the electrons occupy the d xz and d yz orbitals 
creating the antiferro-orbital (alternating orbital) order. 
At the B4 sites only the d xy orbitals are occupied by 
electrons. 



IV. DISCUSSION 

The presented results reveal some interesting aspects 
of two iron-based spinels. As a common feature, the elec- 
tronic properties of both compounds are strongly affected 
by the on-site Hubbard interaction U. The interplay of 
this local Coulomb interaction with a lattice distortion is 
the driving force of the observed phase transitions, which 
lead in both spinels to a new ordered phase with a broken 
cubic symmetry. 

In magnetite, this interaction does not lead directly to 
the insulating phase as in the Mott-Hubbard mechanism. 
It enhances a natural tendency of the degenerated ti g 
states to create an orbital ordering. Such state becomes 



very sensitive to lattice distortion, which by splitting the 
tig band opens the insulating gap. This mechanism is 
similar to the Jahn- Teller effect, however, the lattice dis- 
tortion in magnetite is much more complex than in simple 
perovskites, see e.g. Ref. [28j. The structural transition 
is here induced by two order parameters, which couple to 
each other and generate the monoclinic phase [14J . 

Long time ago, Anderson predicted the existence of 
the short-range order in magnetite on the basis of the 
simple electrostatic model [l|. There are many signs of 
such state above the Verwey transition. One of them is 
the insulating gap that does not close completely at high 
temperatures [29(. Recently, the resonant x-ray scatter- 
ing studies revealed that charge-orbital ordering starts 
to develop about 10 K above Ty [30]. The charge-orbital 
fluctuations couple to phonons producing the critical dif- 
fuse scattering. The neutron studies found the symme- 
tries of lattice distortions associated with these fluctua- 
tions. They are the precursors of the order parameters 
A 5 P and X 3 [13. 

We found that the geologically important mineral 
Fe2SiC>4-spinel is a Mott insulator. The on-site Hubbard 
interaction U in the 3d states on iron generates the in- 
sulating state with the AF order. We suggest that these 
AF interactions drive the phase transition observed at 
Tw = 11.8 K. Furthermore, we discovered that the AF 
order induces the tetragonal distortion reducing the cubic 
Fd3m symmetry to its subgroup Iii/amd. This distor- 
tion releases the geometrical frustration of the magnetic 
interaction enabling for the long-range order. In Ref. 
|19| , we have estimated the exchange constant J and the 
Curie- Weiss temperature #cw using the classical Heisen- 
berg model. The obtained large values of the exchange 
constant J = 1.7 meV and the Curie- Weiss temperature 
&CW — 340 K are in sharp constrast with the low value of 
TV = 11.8 K and indicate that magnetic interactions are 
strongly frustrated and are expected to generate short- 
range order above TV- 

So far, there are no experimental indications that the 
AF order in its simplest form (as in Fig. [SJ exists below 
T/v, and its nature is more subtle. One possibility is that 
the magnetic order is non-colinear, with complicated spin 
orientations. In such a case, the low-symmetry structure 
could be different from the tetragonal. It is also possible 
that below TV the magnetic order has only a short-range 
character (similar to spin ice), or the long-range order 
coexists with the short-range order like in the ZnFe204 
spinel (H). 

In conclusion, we suggest that while the origin of the 
Verwey transition in Fe304 could be understood, the na- 
ture of the magnetic order in Fe2Si04 below Tn poses 
an interesting experimental problem. We hope that fu- 
ture experiments will either confirm the predictions of 
the present studies, or indicate necessary ingredients of 
a more complete theoretical treatment of this problem. 
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